Global regularity for a Birkhoff-Rott-a approximation of the dynamics of vortex 

sheets of the 2D Euler equations 
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I. INTRODUCTION 



One of the novel approaches for subgrid scale model- 
ing is the a-regularizations of the Navier-Stokes equa- 
tions (NSE). The inviscid Euler-a model was originally 
introduced in the Euler-Poincare variational framework 
in In [1,13, SB 01 the corresponding Navier-Stokes- 

a (NS-a) [also known as the viscous Camassa-Holm 
equations or the Lagrangian-averaged Navier-Stokes-a 
(LANS-a)] model, was obtained by introducing the ap- 
propriate viscous term into the Euler-a equations. The 
extensive research of the a-models (see, e.g., [3| . kj, m , |d, 
, 0, [II [13, [II, Q [11 E [II [ll H [2fl, [2ll2l, 
24| ) stems from the successful comparison of their 
steady state solutions to empirical data, for a large range 
of huge Reynolds numbers, for turbulent flows in infinite 
channels and pipes. On the other hand, the a-models can 
also be viewed as numerical regularizations of the origi- 
nal, Euler or Navier-Stokes, systems. The main practical 
question arising is that of the applicability of these reg- 
ularizations to the correct predictions of the underlying 
flow phenomena. 

In this paper we present some results concerning the 
a-regularization of the two-dimensional (2D) Euler equa- 
tions in the context of vortex sheet dynamics. A vor- 
tex sheet is a surface of codimension one (a curve in the 
plane) in inviscid incompressible flow, across which the 
tangential component of the velocity has a jump discon- 
tinuity, while the normal component is continuous. The 
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evolution of the vortex sheet can be described by the 
Birkhoff-Rott (BR) equation [H, [H, HI . This is a non- 
linear singular integro-differential equation, which can 
be obtained formally from the Euler equations assuming 
that the evolution of a vortex sheet retains a curve-like 
structure. However, the initial data problem for the BR 
equation is ill-posed due to the Kelvin-Helmholtz insta- 
bility [25|, i28(] . Numerous results show that an initially 
real analytic vortex sheet can develop a finite time sin- 
gularity in its curvature. This singularity formation was 
studied with asymptotic techniques in 29^, 2Ql Sind nu- 
merically in [30, [3ll [33|. S pecific examples of solutions 
were constructed in [33l.l34*|. where the development, in a 
finite time, of curvature singularity from initially analytic 
data was rigorously proved. 

The problem of the evolution of a vortex sheet can 
also be approached, in the general framework of weak 
solutions (in the distributional sense) of the Euler equa- 
tions, as a problem of evolution of the vorticity, which 
is concentrated as a measure along a surface of codimen- 
sion one. The general problem of existence for mixed- 
sign vortex sheet initial data remains an open question. 
However, in 1991, Delort [s^ proved a global in time 
existence of weak solutions of the 2D incompressible Eu- 
ler equation for the vortex sheet initial data with ini- 
tial vorticity bein g a Radon measure of a distinguished 
sign, see also [36l. ISTI. [ssl . l39l IdOl . j^H. This result was 
later obtained as an inviscid limit of the Navier-Stokes 
regularizations of the Euler equations fSTl . [s^, and as 
a limit of vortex methods [ssl. 140*1 . The Delort's result 
was also extended to the case of mirror-symmetric flows 
with distinguished sign vorticity on each side of the mir- 
ror (4^ . However, the problem of uniqueness of a weak 
solution with a fixed sign vortex sheet initial data is still 
unanswered, numerical evidences of non-uniqueness can 
be found, e.g., in [H, [d^]. Furthermore, the structure of 
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weak solutions given by Delort's theorem is not known, 
while the BirkhofF-Rott equations assume a priori that a 
vortex sheet remains a curve at a later time. A proposed 
criterion for the equivalence of a weak solution of the 
2D Eulcr equations with vorticity being a Radon mea- 
sure supported on a curve, and a weak solution of the 
BirkhofF-Rott equation can be found in ^45*1. Also, an- 
other definition of weak solutions of BirkhofF-Rott equa- 
tion has been proposed in \^ . For a recent survey 
of the subject, see (isj . 

The question of global existence of weak solutions for 
the three-dimensional Euler-a equations is still an open 
problem. On the other hand, the 2D Euler-a equations 
were studied in [i^, where it has been shown that there 
exists a unique global weak solution to the Euler-a equa- 
tions with initial vorticity in the space of Radon measures 
on R^, with a unique Lagrangian flow map describing 
the evolution of particles. In particular, it follows that 
the vorticity, initially supported on a curve, remains sup- 
ported on a curve for all times. 

We present in this paper an analytical study of the 
a-analogue of the Birkhoff-Rott equation, the Birkhoff- 
Rott-a (BR-a) model, which is induced by the 2D Euler- 
a equations. The BR-a model was implemented compu- 
tationally in (soj , where a numerical comparison between 
the BR-a regularization and the existing regulari zing 
methods, such as a vortex blob model 
has been performed. We remark that, unlike the vortex 
blob methods that regularize the singular kernel in the 
BirkhofF-Rott equation, the a-model regularizes instead 
the Euler equations themselves to obtain a smoother ker- 
nel. 

We report in Section HVl our main result, which states 
that the initially smooth self-avoiding 2D vortex sheet, 
evolving under the BR-a equation, remains smooth for 
all times. In this short communication we only report the 
results and sketch some of their proofs, the full details will 
be reported in a forthcoming paper. In Section [IT] we de- 
scribe the BR-a equation. Section [1111 studies the linear 
stability of a flat vortex sheet with uniform vorticity den- 
sity for the 2D BR-a model. The linear stability analysis 
shows that the BR-a regularization controls the growth 
of high wave number perturbations, which is the reason 
for the well-posedness. This is unlike the case for the 
original BR problem that exhibits the Kelvin-Helmholtz 
instability, the main mechanism for its ill-posedness. 



II. BIRKHOFF-ROTT-Q EQUATION 

The incompressible Euler equations in in the vor- 
ticity form are given by 



V ~ K * q, 



(1) 



where K (x) = ^V"*- log |a;I, v is the fluid velocity field, 
q = curl V is the vorticity, and q^" is the given initial vor- 
ticity. 

The 2D Euler-a model [U, i, d, M, M, 113 is an invis- 
cid regularization of the Euler equations, such that the 
vorticity is governed by the system 



-^ + {u-V)q = 0, 
u = K°' * g, 
g(x,0)=(7"X^)- 



(2) 



Here u represents the "filtered" fluid velocity, and a > 
is a length scale parameter, which represents the width 
of the filter. At the limit a = 0, we formally obtain 
the Euler equations H]). The smoothed kernel is X" — 
G" * K, where G" is the Green function associated with 
the Helmholtz operator (/ — a^A), given by 



G" (x) 



a^ Va 



1 1 



Ko 



(3) 



here x = {xi,X2) G and Kq is a modified Bessel func- 
tion of the second kind [58| . 

Let 7W(R^) denote the space of Radon measures on R^; 
Q denote the group of all homeomorphism of R^ , which 
preserve the Lebesgue measure; and rj = ri{-,t) denote 
the Lagrangian flow map induced by ([2|) and obeying the 
equation dtri(x, t) = u(ri(x, t), t), rj {x, 0) = x. 

Oliver and ShkoUer H9\ showed global well-posedness 
of the Euler-a equations ^ with initial vorticity in 
A4(R^) (which includes point- vortex data). 

Theorem 1. (Oliver and ShkoUer ^^\) For initial data 
q™ G A^(R^), there exists a unique global weak solution 
(in the sense of distribution) to ^ with 

V e G^ (R;a), u e G(R;G(R2)) , q e C {R; M{m.^)) . 

The Birkhoff-Rott-a equation, based on the Euler-a 
equations is derived similarly to the derivation of the orig- 
inal Birkhoff-Rott equation. Detailed descriptions of the 
BirkhofF-Rott equation as a model for the evolution of the 
vortex sheet can be found, e.g., in [S^, |4l|, [5§|. We re- 
mark, that while the BR equations assume a priori that 
a vortex sheet remains a curve at a later time, in the 2D 
Euler-a case, if the vorticity is initially supported on a 
curve, then due to the existence of the unique Lagrangian 
flow map given by Theorem [1] it remains supported on 
a curve for all times. Hence the BR-a equation gives an 
equivalent description of the vortex sheet evolution, as 
the 2D Euler-a equations. It is described in the follow- 
ing proposition. 

Proposition 2. Let q be the solution of ([2]) in the 

sense of the Theorem [7[ Assume, furthermore, that q 
has the density ^{a,t) supported on the sheet (curve) 
j:(t) = {x = x{a,t) eR'^\ao{t) <a <ai{t)}, that is, 
the vorticity q{x,t) satisfies 

ip{x)dq{x,t) ^ / ip{x{a,t))-f{a,t)\xa{<J,t)\da, 
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for every ^ € (R^), 7(-,t) G L'^{\x^\da). Then this 
sheet evolves according to the equation 



-Q-x{a.t) 

<Tl(t) 



if" (x {a, t) - x{a',t))-f (ct, t) \xa (cr', t)\da' . 



Additionally, if T (cr, t) — J^_, 7 (cr', t) \xa (cr', i)| dcr', 
where x (cr*,t) is some fixed reference point on S(t),, de- 
fines a strictly increasing function of a (e.g., as in the 
case of positive vorticity), then the evolution equation is 
given by the Birkhoff-Rott-a (BR-a) equation 



d_ 

dt 



x{T,t) 



ri 



K°'(x{V,t)-x{T',t))dT' (4) 



with 7 = l/|2^r| being the vorticity density along the sheet 
and —00 < Fq < Fi < 00. 

Here cro,ai can represent either a finite length curve, 
or an infinite one. Notice that 



where 



if" (x) = V^*" (|a;|) = — tD*" (bl) , 

a; 



(r) = — 



and 



^0 ( 


— j + logr 










1 












2^ 







D^<^{r) = -— (r) 
dr 



Kq and Ki denote modified Bessel functions of the sec- 
ond kind of orders zero and one, res pec tively. For de- 
tails on Bessel functions, see, e.g., 58|. We remark 
that the smoothed kernel (x) is a bounded contin- 
uous function, that for ^ behaves as K°' (x) — 

— ^-^x-^\os— + O ( ^-^] . That is, it is non-singular 

kernel. Since j{-,t) G L^{\xa\da) we can show the in- 
tegrability of the relevant terms, even though jif" {x) \ is 
decaying like "'^ at infinity. 



III. LINEAR STABILITY OF A FLAT VORTEX 
SHEET WITH UNIFORM VORTICITY DENSITY 
FOR 2D BR-a MODEL 

The initial data problem for the BR equation is highly 
unstable due to an ill-posed response to small perturba- 
tions called Kelvin-Helmholtz instability [2^, |28|. The 
linear stability analysis of the BR-a equation shows that 
the ill-posedness of the original problem is mollified, and 
the Kelvin-Helmholtz instability of the original system 
now disappears. 



When the vortex sheet can be parameterized as a graph 
of a function in the form X2 — X2 {xi , t) the BR-a system 
(|4]) takes the form 



dx2 
l)t 



dx2 

"1 + U2, 

0x1 



(5) 



with velocity u = (ui, U2) given by 
u (xi, t) = p.v. / K"" {x{xi,t) ~ x{x'i,t))j{x[,t)dx'i, 



where x{xi,t) = {xi,X2 (a;i,i))*. The fiat sheet X2 = 
with uniformly concentrated intensity 70 is stationary so- 
lution of ([5]). By linearization about the fiat sheet we 
obtain the following linear system 



dx2 



U2, 



dui 



dt ^'^'dxi 



where 



ui {xi,t) = -70 {sgn{xi)D'^°' {\xi\)) * 
U2 {xi,t) = (sgn(a;i)L>*"(|a;i|))*7, 



dX2 

dx\ ' 



and (i2,7) is a small perturbation about the fiat sheet. 

Consequently, the equation for the Fourier modes is 
given by 



d_fx2 
dt U 

where 







-iW sgn(fc)d(fc) 



:Sgn(A:)d(fc)\ /^^ 



(6) 



d{k) 



1 



1 



-1/2 



1. 



Observe that in order to calculate the Fourier transform 



(sgn (a;i) Z?*" (|a,i |)) {k) = - sgn(fc)d(fc), 

we used the integral representation of the modi- 
fied Bessel function of the second kind Ki {xi) = 

xi/^°°e-=^i* (t2 - (see, e.g., H). The eigen- 

values of the coefficient matrix, given in are 



A(fc) = ±-|7o||fc| 1 



1 



x^k"^ 



-1/2N 



(7) 



To conclude, the a-regularization mollifies the Kelvin- 
Helmholtz instability as follows: we have an algebraic 
decay of the eigenvalues to zero of order ^jsj^ , as /c ^ 00 
(a fixed). While, for a ^ 0, for fixed fc, we recover the 
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eigenvalues of the original BR equations ±i |7o| |fc| (see, 
e.g., M)- 

For the sake of comparison, we observe that for the 
vortex blob regularization of Krasny [32| , where the sin- 
gular BR kernel, K{x), was replaced with the smoothed 
kernel 

|2 



Ks (x) = K (x) 



1 



(52 27r 



the eigenvalues are 



A(fc) = ±-e 



-Sk 



7o||fc| 



with an exponential decay to zero, as fc — > oo ((5 > 
is fixed). As (5 ^ 0, for fixed k, one recovers again the 
eigenvalues of the original BR equations. 

The behavior of the eigenvalues of the linearized sys- 
tem ^ indicates that high wave number perturbations 
grow exponentially in time with a rate that decays to 
zero, as fc — !• oo, which is the reason for well-posedness 
of the a-regularized model. This is unlike the original 
BR problem that exhibits the Kelvin-Helmholtz insta- 
bility. It is worth mentioning that the a-regularization 
is "closer" to the original system than the vortex-blob 
method at the high wave numbers, due to the alge- 
braic decay instead of exponential one in the vortex blob 
method. This result was also evaluated computationally 
in [H. 



IV. GLOBAL REGULARITY FOR BR-q 
EQUATION 

In this section we present the global existence and 
uniqueness of solutions of the BR-a equation ^ in the 
appropriate space of functions. We show that initially 
smooth solutions of ^ remain smooth for all times. 

Let us first describe the Holder space 
C"''3 (E C K; ^ < /5 < 1, which is the space 
of functions x : S C M ^ K.^ . with finite norm 



C"./3(E) 



E 

*;=0 



CO(E) 



/3(S) 



where 



lco(E) = sup |a;(r)| 
res 



and |-|^ is the Holder semi- norm 



\x{T)-x{T')\ 



l/3(E) 



sup 

r.r'es |r — T' 



We also use the notation 

= inf 
rj'GS 



\x{r)-x{r')\ 
ir-r'i 



Next we state our main result. 



Theorem 3. Let n > 1, < f3 < 1, 

x(r,0) = xo{T) E C^'^{To,Ti) n {\x\^>0}, 
then for any T > there is a unique solution 
xeC'{[-T,T];C^'0{To,T^)n{\x\^>O}) of Q. 
In particular, if xq G C°° (ro,ri) n {|a;|^ > 0} then 

xeC ([-r,r]; (ro,ri) n {\x\, > o}). 

We remark that, although the kernel K" is a continu- 
ous bounded function, its derivatives are unbounded near 
the origin, and the condition |a;|^ > 0, which generally 
means self-avoiding curves, allows us to show the integra- 
bility of the relevant terms. Furthermore, it is also worth 
mentioning that being bounded away from zero is 
similar to the chord arc hypothesis [6l|, used later in 
Hill . 

Now we sketch the main steps involved in the proof 
of Theorem [3l First, we apply the Contraction Mapping 
Principle to the BR-a equation ^ to prove the short 
time existence and uniqueness of solutions in the appro- 
priate space of functions. We show that initially C^'^ 
smooth solutions of (|4]) remain C^''^ smooth for a finite 
short time. Next, we derive an a priori bound for the 
controlling quantity for continuing the solution for all 
time. Then we extend the result for higher derivatives. 
The full details will be reported in a forthcoming paper. 

Sketch of the proof. We consider the BR-a equation as 
an evolution functional equation in the Banach space 

^ (F, t) = ll^^ [x (F, t) X (F', t)) rfF', 
x(F,0) = 2;o (F) eC"''^n{|a;|, >0} 

with 7 — l/|a;r| being the vorticity density along the 
sheet. Notice that the initial density is well defined for 
the subset {\x\^ > 0}. 

Step 1. We show the local existence and uniqueness of 
solutions. To apply the Contraction Mapping Principle 
to the BR-a equation ^ we first prove the following 
proposition 

Proposition 4. Let 1 < M < oo, — oo < Fq < Fi, and 

let 5*^ be the set 

|f K-> X (F) e (Fo, Fi) , \xr\co < M, \xl > ^| . 
Then the mapping x (F, t) i— > 



u{x{T,t),t)^ / K°' {x{T,t) - x{r',t))dT' 

defines a locally Lipschitz continuous map from S^^ into 

This implies the local existence and uniqueness of so- 
lutions: 



5 



Proposition 5. Given x^{T) G C^^f^ {Tq,Vi) D 
{|a;|^ > 0}, there exists 1 < A/ < oo and a time T{M), 
such that the system ([8]) has a unique local solution 
X e Ci((-T(M),T(M));5*0- 



Step 2. The obtained local solutions can be continued 
in time provided that we have global, in time, bounds 



and |xr(-,t)|^. To control these quantities 



dt. We 



we need to bound \\\/xu{x{-,t),t)\ 
sketch the proof of this bound. We write WxU {x(r, t), t) 



as 



Vxuix{T,t),t) 



V,if" {x {T,t)~x{T',t))dr' 



= h+h, 



where 



i;.Hr'c(r„,r.): l-^"--"-g''"l <. 



for a fixed small < e < 1, to be further refined later. 
Let r] denote the unique Lagrangian flow map given by 
Theorem [TJ Denote the distance between two points 
r]{x,t) and r](x',t) by r (t) = |ry (a;, i) — ry (a;', where 
r(0) = |a:-a;'|. 

Then, using the estimate (2.14) of [4^, we have 



where Ci 



^ lk"llA^(K2)- Now, using also that \xo\ 



is bounded away from zero, we can bound 
from below, which in turn implies the bound 

/i <c(i,^,||'Z™||^(R.),|a;oL) . 

While to bound I2, we use the boundness of 
|V,if"(x(r,i),a:(r',i))| in {F' e (rcTi) : 

\x{r,t)-x(^r' ,t)\ 



> s}. Hence 



\Vxu{x{-,t),t)\\ 



i°°(ro,ri) 



dt < 



<C\-,T,\\r\\j^^^.y\xo\, 



(10) 



Now, by Gronwall inequality the bound (jlOp provides 
and \xr {■,t)\^o on [0, T]. The bound 



bounds on j^rjji ^nu |xr v,oj\co 
on |a;r {■,t)\p on [0,T] is a consequence of 



dt 



xr (r, t) = {x (r, t) , i) • xr {T, t) , 



|V,u(a;(-,t),i)|^ < c(^i,|xrL=o ,|xL,ri-ro^ 



pO|l and the Gronwall inequality. 

This yields global in time existence and uniqueness of 
C^'P solutions of m. 



dt 



where 



r(t) < / \K''{x,y)-K''{x',y)\\q{y,t)\dy 

-cU C-^) \W-\\ 



0, r = 0, 
(p[r) — r (1 — logr) , < r < 1, 

1, r>l. 



Step 3. To provide an a priori bound for higher deriva- 
tives in terms of lower ones, we show that for x G 
gM ^(jn,f3 (ro,ri), 

I" {x {■,t) , 0|„,/3 < C (^i, M, \x (., t)l_^ J^ \x (., t)l^p , 

hence by Gronwall inequality and the induction argu- 
ment, it is enough to control |a^(-,i)|^ and |a::r(-,i)|^, 
to guarantee that x{T,t) £ C"^'' (Fo, Ti), for all n > 
1, (and consequently in C°°(ro,ri), whenever xq G 

c-(ro,ri)n{|a;L >o}). □ 



By comparison with the solution of the differential equa- 
tion 



r{t) = ~C-^ ) iir 

at a \ a 



Im(r2) 



we can choose s 



F (t - llo*"ll 

, , r \x(r,t)-x(r' ,t)\ 
such that, for — - — !^ < e, 



small enough, 



\x{T,t)~x{T',t)\ 



> 



a 



> 



|x(r,o)-x(r',o)| 



(9) 



V. CONCLUSIONS 

The 2D Euler-a model P, i, H, HE [H, [13 is an invis- 
cid regularization of the Euler equations. In [3^ it has 

been shown the existence of a unique global weak solu- 
tion of 2D Euler-a equations, when the initial vorticity 
is in the space of Radon measures on R^. The Birkhoff- 
Rott-a equation for the evolution of the 2D vortex sheet 
is induced by the 2D Euler-a equations, and it is an a- 
analogue of the Birkhoff-Rott equation, induced by the 
2D Euler equations. 

The structure of weak solutions of 2D Euler equations, 
for the vortex sheet initial data with initial vorticity being 
a Radon measure of a distinguished sign, given by Delort 
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di, m, [13, H Hi, \M, mi is not known, yet the BR 
equations assume a priori that a vortex sheet remains a 
curve at a later time. On the contrary, in the 2D Euler-a 
case, if the vorticity is initially supported on a curve, it 
remains supported on a curve for all times, hence the BR- 
a equation gives an equivalent description of the vortex 
sheet evolution, as the 2D Euler-a equations. 

In this paper we report the global regularity of the 
BR-a approximation for the 2D vortex sheet evolution. 
We show that initially smooth self-avoiding vortex sheet 
remains smooth for all times, under the condition that 
the initial density is an integrablc function of the vortex 
curve with respect to the arc-length measure. 

Unlike the original BR problem that exhibits the 
Kelvin-Helmholtz instability, the linearized, about the 
flat solution, BR-a model has growth rates that decay 
to zero for large wave numbers, larger than 0{a). This, 
in turn, is also an indication of the role that the param- 
eter a plays in slowing the process of formation of scales 
smaller than a. Another indication that a controls the 
development of small scales, smaller than a, arises from 
the Lagrangian description of the flow. The lower bound 
^ implies that the evolution of small scales, relative to 



a, at each instant of time, is controlled from below by 
the initial ratio. That is, for any finite time, the spatial 
scales smaller than alpha develop at a controlled rate. 

The linear stability analysis also implicates that the 
BR-a approximation could be closer to the original BR 
equation than the existing regularizing methods, such as 
vortex blob model, due to the less regular kernel. A nu- 
merical study comparing the a and the vortex blob regu- 
larizations for planar and axisymmctric vortex filaments 
and sheets is reported in 50]. 

The full details of the results reported in this paper 
will be presented in a forthcoming paper. 
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